1-D Mott insulator transition of a Bose-Einstein condensate 
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The superfluid to one-dimensional Mott-insulator transition of a ^^Rb Bose-Einstein condensate 
is demonstrated. In the experiment, we apply a one-dimensional optical lattice, formed by two 
laser beams with a wavelength of 852 nm, to a three dimensional BEG in a shallow trap. We use 
Kapitza-Dirac scattering to determine the depth of the optical lattice without knowledge of its exact 
geometry. We further study the dynamics of the transition as well as steady-state phase behavior 
specific to the one- dimensional case. 
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Bose-Einstein condensation (BEG) was first demon- 
strated in 1995 [1, 2], creating an explosion of interest 
in previously unattainable many-body quantum phenom- 
ena. Of particular interest is the phase transition from 
superfluid to Mott insulator [3]. In the general, three- 
dimensional (3-D) case, this phase transition occurs when 
a 3-D optical lattice is applied to a BEG. As the lattice 
depth is increased, the BEG transitions from a super- 
fluid state to a state with a definite number of atoms 
in each lattice well. The Mott insulator transition has 
drawn the interest of both atomic and condensed matter 
physicists, due to the possibilities it creates for simulat- 
ing ideal, controlled condensed matter systems. Under 
the correct circumstances the transition could be used 
to create supersolids or other novel phases of matter j^. 
Doping the Mott insulator with fermions can be used 
to simulate a semiconductor [5]. One can use Feshbach 
resonances to create molecules in a Mott insulator with 
two atoms per site l6j, which could eventually lead to a 
molecular BEG 0, 0. The Mott insulator state could 
also provide a means to entangle neutral atoms and form 
a quantum register for a quantum computer (for a re- 
view, see Ref. [9]). Several laboratories have succeeded 
in producing the transition from a BEG to a 3-D Mott 
insulator [3, 5, 6, 10]. In lower dimensions, the Mott in- 
sulator transition has been achieved using; 2-D p m 
and 1-D p^] Bose gases. One of the difficulties in at- 
taining the Mott insulator transition is that the lattice 
depth necessary to induce the transition requires high 
laser power or, more commonly, narrowly focused beams. 
Demonstrating the transition from a BEG to a 1-D Mott 
insulator (a 3-D BEG in a 1-D lattice formed by two laser 
beams) is even more difficult because a deeper lattice is 
needed jl^^llll- Furthermore, it can be difficult to deter- 
mine the exact lattice depth due to uncertainties in the 
lattice beam alignment. In this paper, we experimentally 
investigate the properties of Mott insulator transition of 
a 3-D BEG in a 1-D lattice using a novel, robust method 
to calibrate the lattice depth. In contrast to Ref. [HI , we 
apply no transverse lattice potential [16]. The quantum 
gas therefore retains superfluidity in the two dimensions 
transverse to the lattice-beam direction, with observable 
consequences. 

In the experiment, we start with a ^^Rb BEG of 5- 
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FIG. 1: (a) TOF images of Kapitza-Dirac scattering as a 
function of lattice duration At, in steps of 0.5 /is, starting 
from 0.5 /xs. (b)-(d) Scattering ratios Nu/Nt, as defined in 
the text, for the scattering orders n = 0,1, and 2, respectively, 
as a function of lattice duration. The data in (b) are fit with 
Jo (a At), with best-fit parameter a=0.79. The lines in (b)-(d) 
show Jn(0.79At) with respective values of n. 



8x10^ atoms in a practically harmonic magnetic trap. 
Our BEG apparatus is described in [T3]. In the data 
presented, we use magnetic-trap frequencies of 80 Hz and 
200 Hz in the direction of the optical lattice (in the other 
two directions, the trap frequencies are 20 Hz and 80 Hz 
or 40 Hz and 200 Hz, respectively). The optical lattice is 
formed by a retro-reflected, far-off-resonance laser beam 
(wavelength 852 nm, power up to 200 mW after flber). 
The beam is focused into a spot with an intensity full- 
width half-maximum of 80 /im. The depth of the optical 
lattice is obtained as follows. 

When a standing wave of sufficiently short duration is 
applied to cold atoms such that the atoms are stationary 
while the lattice is applied, the system is in the Kapitza- 
Dirac scattering regime (analogous to the Raman-Nath 
regime in optics) . The 1-D optical lattice adds a potential 



V(x) = -Fo(1-cos(2A:lx)) 



(1) 



to the atoms over the time interval At that the lattice is 
on. 2Vo is the lattice depth, and /cl is the wavenumber 
of the lattice beams. Assuming an initial wavefunction 
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ip{x,t = 0) = 1, the wavefunction after the lattice pulse 
is, neglecting a global phase factor, 



ip{x^ t > 0) = exp ( i 
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W"^n f ^jexp(z2nfcL:r) . 



(2) 



The expression in the sum shows that the BEC breaks up 
into momentum components that are integer multiples 
of 2hki, with amplitudes given by Bessel functions. In 
particular, the order n = first vanishes at a time Ato for 
which Vo = 2A048h/Ato. The lattice depth 2Vb can thus 
be found by measuring the time At at which the n = 
order first vanishes. For an atomic polar izability a and a 
single-beam lattice intensity /i, the lattice depth is also 
given by 2Vb = a/i/(2ceo). Using the above equations, 
the lattice depth can be experimentally calibrated against 
arbitrary linear functions of /i, such as the measured 
beam power. The strength of this method is that no 
geometrical measurements of the lattice beam size and 
position are necessary. 

To examine the Kapitza-Dirac scattering, the optical 
lattice is applied to the BEC for a few microseconds. The 
lattice and the trap are then turned off simultaneously. 
The BEC is allowed to expand freely for times of flight 
(TOP) of 16 ms or 12 ms, for the case of the 80 Hz 
or 200 Hz traps, respectively. After the expansion, we 
take absorption images, shown in Fig. [1] (a). Since the 
BEC temperature depends on the trap frequency, differ- 
ent trap frequencies require different TOFs to produce 
the best image. For each scattering order n, we measure 
the atom number {Nt is the total atom number). 
In Fig. [U (b) we plot Nq/Nt vs At and find the lattice 
duration Ato where this ratio first approaches zero; the 
corresponding lattice depth is 4.8096fi/Ato. In this way, 
we can determine the proportionality constant relating 
the lattice-power reading to the lattice depth. Typically, 
we find Ato ~ 3 /is, which is well within the validity 
range of the short-pulse approximation underlying the 
above treatment, as evidenced by the fits in Fig. [1] (b) 
which match the data well up to 6 /is. In the measure- 
ments presented this paper, this calibration procedure is 
repeated frequently to account for small changes in the 
alignment of the lattice beams. Furthermore, this pro- 
cedure is used to initially align the lattice with respect 
to the BEC: better overlap between the lattice and the 
BEC leads to a smaller Ato- 

With the lattice depth calibrated, we investigate the 
transition from a superfluid to a 1-D Mott insulator. 
While in Kapitza-Dirac scattering a sudden, microsec- 
ond pulse is applied, the transition of the BEC from a 
superfluid to a Mott insulator requires adiabatic load- 
ing of the atoms into the optical lattice, taking of order 
ten milliseconds. We use the amplitude modulation of 
an AOM to control the power of our lattice beam. We 
ramp up the power of the lattice beam from zero to its 
final value over 10 ms, and then hold it there for 5 ms. 
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FIG. 2: Left: TOF images. Right: lattice depth as a function 
of time for (a) BEC with no lattice, (b) superfluid phase, (c) 
1-D Mott insulator, (d) superfluid phase recovered after Mott 
insulator, and (e) BEC with no lattice, recovered after Mott 
insulator. 



Then, we simultaneously turn off the lattice and the mag- 
netic trap, and take a TOF measurement. As can be 
seen in Fig. [2] (b), for small lattice depths the BEC is 
only slightly modulated by the lattice, corresponding to 
the appearance of only two weak side peaks, at ±2fi/cL- 
As the depth of the lattice is increased, the transition 
to Mott insulator occurs. The typical signature of the 
transition is that the side peaks disappear and the cen- 
tral peak broadens, reflecting the momentum distribution 
of the localized wavefunction in a single lattice well 
Here, we find that the system fully reaches the 1-D Mott 
insulator state around 30 £^rec- 

The Mott insulator transition is a quantum phase tran- 
sition, and thus is reversible; to be certain that we have 
seen the Mott transition, as opposed to a lattice-induced 
heating effect, we must show that we can reverse it. To 
demonstrate this, we ramp the lattice to 31 Erec over 
10 ms, hold it there for 5 ms, and then ramp back down 
over 10 ms. As can be seen in Figs. [2] (d) and (e), we 
obtain a modulated superfluid and BEC when we ramp 
down to a weak lattice and no lattice, respectively. Thus, 
the effect we see is fully reversible, providing strong evi- 
dence that it is the 1-D Mott insulator transition. 

In the 1-D Mott-insulating state (Fig.[2](c)), the quan- 
tum gas loses phase coherence in the direction of the op- 
tical lattice while retaining its superfluidity in the other 
two directions. The 1-D Mott insulator can thus be 
thought of as a stack of uncorrelated pancake BECs, 
each containing ~3000 atoms under the conditions of 
Fig. [2] (c). As can be seen in Fig. [2] (c), the 1-D Mott 
insulator expands much farther in the direction of in- 
sulation than in the directions of superfluidity. This 
is largely due to the momentum spread of the pan- 
cake BECs in the lattice-beam direction. Examining the 
TOF image in Fig. [2] (c), we find a velocity spread of 
Ap/mRb = 8 mm/s. Using the Heisenberg uncertainty 
relation, AxAp > this corresponds to a localization 
Ax = 46 nm, or 11% of the lattice period. Neglect- 
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ing mean-field effects and using tiie fact tiiat tiie lattice 
wells are approximately harmonic near their minima, we 
find an oscillation frequency of 27r x 35 kHz for a lattice 
with a depth of 30 Ej-ec^ and velocity and position un- 
certainties of 8.9 mm/s and 41 nm, respectively, for the 
ground state. These numbers match the values derived 
from Fig. [2] (c) quite well, showing that the expansion in 
the lattice-beam direction is mostly driven by the kinetic 
energy of the pancake BECs in the optical-lattice wells. 

A more subtle effect is that in the insulating case 
the expansion transverse to the lattice-beam direction is 
considerably slower than in the lattice- free BEC: about 
1.5 mm/s and 2.5 mm/s, respectively. We attribute the 
difference to a variation in the manifestation of the re- 
pulsive mean-field potential (estimated to be < 1 kHz 
for our BECs in 200 Hz magnetic traps). Without the 
lattice, the BEC expansion is driven by a combination 
of the mean-field pressure and the kinetic energy of the 
BEC in the magnetic trap, leading to a final expansion 
speed of about 2.5 mm/s in all directions in Fig. [2] (a). 
After application of the deep, Mott-insulating lattice in 
Fig.[2](c), the expansion is mostly driven by the compar- 
atively high kinetic energy of the BEC pancakes in the 
optical-lattice wells, leading to a much faster expansion 
in the lattice direction. The faster expansion leads to 
a reduction of the time over which a substantial mean- 
field pressure exists, leading to a reduced final expansion 
speed transverse to the lattice, as observed. 

To quantitatively characterize the 1-D Mott transition, 
we examine it as a function of lattice depth. We use the 
visibility of the side peaks, to map out the transition: 



Na-Nb 
Na^Nb 



(3) 



where A^a is the linear atom density of one side peak, 
and A'b is the linear atom density at the minimum be- 
tween the center peak and the side peak. The timing of 
the lattice application is as described above. Examining 
the resulting image, we take the linear atom density as 
a function of position in the lattice direction along the 
central strip of the image, integrating over three pixels 
in height (a pixel in the image corresponds to 6.7 /im). 
For Na^ we choose the local maximum at the side peak, 
if there is one, and for Nb we choose the local minimum 
between the side peak and the central peak, as shown 
in the sample data in Fig. [3] (a). If there is no local 
side maximum, we designate v = 0. We calculate the 
visibility separately for the side peaks on the left and 
right, and repeat the calculation for five separate images 
at each lattice depth. We then average the ten resulting 
values to get the visibility plotted in Fig. [3] (b). As can 
be seen, the Mott transition starts around 10 £^rec7 where 
V first reaches a value lower than its initial value. The 
BEC has fully transitioned to the Mott insulator state 
by 30 ^rec, where v = 0. Our system contrasts strongly 
with that of a 1-D Bose ^as, where the Mott transition 
is complete around 10 ^rec UM- 3-D, the Mott tran- 
sition is complete around 20 ^rec [3|. This is in general 
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FIG. 3: (a) Linear atom density distribution for a BEC after 
12 ms TOF, released from a lattice with a depth of 9.3 ^rec 
and magnetic trap with 200 Hz frequency. Na is the height of 
the side peak and Nb is the height of the valley, (b) Visibility 
as a function of lattice depth from a 80 Hz magnetic trap 
(circles) and a 200 Hz trap (triangles), (c) TOF images as a 
function of lattice depth for a 200 Hz magnetic trap. 



accordance with the prediction in Ref. [14] that the 1-D 
Mott transition requires a deeper lattice than the 3-D 
Mott transition. The specifics do not match, however, in 
that the lattice depth predicted in Ref. [13] necessary for 
the Mott transition is upwards of 50 ^rec for a system 
like ours with ~3,000 atoms per lattice well. As can be 
seen in Fig.[3](b), the Mott transition happens under the 
same lattice conditions in the 80 Hz trap and the 200 Hz 
trap. This indicates that at both 80 Hz and 200 Hz the 
trap has no noticeable effect compared to the lattice. 

Ideally, the visibility should be unity until the Mott 
transition starts [18.]. In the experiment, however, even a 
minute thermal background will lower the visibility, and 
will disproportionately affect images with lower Na- As 
Na^ Nb — ^ 0, this will prevent the visibility in Eq. [3] 
from reaching unity. For these reasons, the first few data 
points in Fig. [3] (b) not only fail to approach unity, but 
are even lower than visibilities seen for lattice depths ~ 
5 E^Qc- We are, however, confident that these issues do 
not affect our determination of where the Mott insulator 
transition occurs, because the transition happens at a 
lattice depth where both Na and Nb are large. 

When the lattice is ramped up, the pancake BECs in 
the individual lattice wells maintain a global phase until 
the Mott insulator transition is reached. Even beyond 
that point, the pancake BECs in their separate wells re- 
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FIG. 4: (a) TOF images as a function of holding time in 
a lattice of 23 Erec depth in a 200 Hz magnetic trap, (b) 
Visibility as a function of holding time of the lattice, for an 
80 Hz magnetic trap (circles) and a 200 Hz trap (triangles). 




FIG. 5: TOF images of Mott insulators from magnetic traps 
with frequencies (a) 80 Hz and (b) 200 Hz. The speckles vary 
from shot to shot, and are due to interference during TOF of 
pancake BECs from different lattice wells. 



quire a certain time to entirely dephase with respect to 
each other; only after the dephasing can one see the typi- 
cal signatures of the Mott insulator transition (for exam- 
ple, in Fig. [3] (c) for lattices of depths larger than about 
20 E^Qc)- For deeper lattices, the time needed to dephase 
is too short to observe. Within the range of the superfluid 
to 1-D Mott insulator transition, however, the dephasing 
time can become sufficiently long to be observable. To 
measure the dephasing time, we start by ramping the 
lattice over 10 ms to 23 £^rec5 leave the lattice on for a 
variable hold time, and take TOF images as a function 
of the hold time, as shown in Fig. [4] (a). The visibilities 
in Fig. [11(b) are obtained as described above. In the case 
of Fig.m the dephasing takes of order three milliseconds. 
In many of the images of the Mott insulator, bright 



and dark vertical stripes appear. These stripes can be 
seen in nearly all of the images of the Mott insulator 
from the 80 Hz trap, and several, but not all, of the 
Mott insulator images from the 200 Hz trap. For exam- 
ple, in Fig. [3] (c), the stripes are visible in the images 
at 23 and 35 £^rec- The arrangement of the stripes ap- 
pears random, with no repetition of the pattern from 
image to image, but the characteristic size of the stripes 
remains the same for a given magnetic trap frequency. 
We believe that these stripes represent interference be- 
tween pancake BECs from different lattice wells during 
TOF. In the 1-D Mott insulator state, each pancake BEG 
still has a definite, but random phase. During TOF, the 
pancakes all expand into each other, leading to a char- 
acteristic interference speckle size. In Fig. [5] we show 
interference patterns in TOF images for the two different 
magnetic trap depths. To find the characteristic speckle 
size. As, in these images, we take fast Fourier transforms 
(EFT) of five images and average them for each trap 
depth. The value of As is given by the inverse spatial 
frequency where the EFT signal reaches the noise fioor. 
We find As = 17 /im for the 80 Hz trap and 27 jam for 
the 200 Hz trap. Using a straightforward analysis, we 
estimate that the number of interfering pancakes, P, is 
related to As and the TOF, T, via P = 2hT/{mAsX). 
From this we find P = 20 and 10 pancakes for the 80 Hz 
and 200 Hz traps, respectively. These numbers agree 
reasonably well with our measurements of the size of the 
BEC in the respective magnetic traps. 

In conclusion, we have demonstrated the 1-D Mott in- 
sulator transition by applying a 1-D optical lattice to a 
BEC. We used Kapitza-Dirac scattering as an accurate 
way to calibrate our lattice depth. We examined time-of- 
fiight images of the BEC as a function of lattice depth, 
and found that the BEC is fully in the Mott insulator 
state at 30 £^rec, and requires a certain time to dephase. 
Furthermore, we observe random interference patterns 
in time-of-fiight images, providing evidence that the 1-D 
Mott insulator consists of pancake BECs that still retain 
superfiuid properties, including a definite phase. 
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